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We suggest measuring the one-particle density matrix of a trapped ultracold atomic cloud by 
scattering fast atoms in a pure momentum state off the cloud. The lowest-order probability for 
the process, resulting in a pair of outcoming fast atoms for each incoming one, as well as of its 
time reversed counterpart, turns out to be given by the Fourier transform of the density matrix. 
Accordingly, important information about quantum correlations can be deduced directly from the 
differential scattering cross-section of these processes. Several most interesting cases of scattering 
- from a single condensate containing a vortex, and from a split condensate characterized by some 
phase difference - are discussed. 
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Achievements of collective quantum states in confined clouds of alkaline atoms and in atomic hydrogen make 
possible studying quantum coherent properties of these systems as well as the revealing fundamental kinetic processes 
leading to the formation of the coherence. The coherence of the condensate has been demonstrated in the atomic 
interference fringes j|. Recently Bragg spectroscopy has been successfully employed to establish that the coherence 
length of a confined condensate is of the order of the cloud size [|J . 

The mechanism of formation of quantum correlations is a matter of great attention and controversy. The emergence 
of the condensate is associated with the formation of off-diagonal long-range order (ODLRO) ||. A primary object 
displaying such an order is the one-particle density matrix (OPDM) p(xi,X2). Typical distances over which these 
correlations become important are comparable with the interatomic spacing r a . Consequently, an "early" detection 
of such emerging correlations is very difficult to achieve by light with the wavelength A>r a usually employed for 
, probing the cloud density. Information about short-range density correlations (at distances < r a ) can, in principle, 

■ be obtained from the absorption of detuned resonant light || . The change in local m-body density correlations (the 
C^l , so-called m!-effect 0) can be seen by measuring recombination rates, and this has been already done experimentally 

for the equilibrium case ||. However, measurements of the correlation length r c (which is far less than A) of the 
forming ODLRO seem very unlikely to be achieved by these methods. Thus it is tempting to have a tool which could 

■ make possible seeing the OPDM directly without fundamental limitations on the accessible distances. 

Q\ , In this paper, we suggest a method of detecting the OPDM which relies on the scattering of coherent atoms off the 
atomic cloud. 

The methods of scattering neutrons and He atoms ]TT| off liquid He are well known. For an incident neutron 

(or atom) which has a relatively small kinetic energy if compared with a typical potential energy of atoms in the 
target, a typical scattering event results in the creation of collective excitations. In contrast, an incident fast neutron 
(atom) ejects an almost free fast atom from the liquid. While the former scenario has practically the same features 
in both helium and weakly interacting Bose gases, the latter has very different consequences in the two systems. 
O ■ Indeed, the weakly interacting Bose gas is almost transparent for the fast atoms, provided the cloud size R traversed 
\ by them is less than the free-path length I = 1/na 2 , where n and a stand for the density and the scattering length, 
respectively. We, however, note that this condition is by no means fundamentally restrictive as long as / is much 
larger than the interparticle distance n" 1 / 3 . In fact, the condition I ^s> n -1 / 3 is guaranteed by the smallness of the 

■ gas parameter £ = na 3 which can be as small as £ ~ 10~ 5 in the atomic clouds. Thus, in a weakly interacting gas 
there is the opportunity of measuring directly the outcome of a single scattering event. We demonstrate that such a 
measurement, if accurate enough, yields direct information on the OPDM in the Bose gas. This method can also be 
applied to studying trapped Fermi gases. In contrast, this is impossible in principle in liquid helium where I « n -1 / 3 
which leads to dominance of the multiparticle excitations in the final-state channel. Accordingly, only some integral 
of the population factors can be obtained from the scattering cross-section || . 

Here, we consider the fast incoming atoms to be identical to the atoms of the cloud. A fast atom carrying momentum 
k can transfer a substantial part of its energy to one atom in the cloud, so that in the final state there are two fast 
outcoming atoms with momenta ki and k2. The quantum-mechanical probability of this process in the lowest order 
turns out to be proportional to a Fourier transform of the OPDM. Let us call this process direct. This process has been 
considered in Ref. \L^ . There is another process which is connected by time reversal with the direct one. Specifically, 
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two fast incoming atoms carrying the respective momenta ki and k 2 enter the cloud and experience stimulated (by 
the atoms of the cloud) scattering, so that in the final state there is only one fast atom with the momentum k. We 
will call this process reversed. 

It should be mentioned that there is another process of atomic scattering, which is analogous to the elastic scattering 
of light. Its probability depends on the cloud density insensitive to r c |l5| ]. 

An expression for the cross-section of the direct and the reversed processes can be derived in the lowest order with 
respect to the two-body interaction for the simplest case of spin-polarized bosons. The essential requirement is that 
the kinetic energy k 2 /2m (we employ atomic units so that Ti — 1) of a fast atom of mass m is much larger than a 
typical interaction energy per particle of the densest p art of t he cloud. This typical energy is given by the chemical 
potential fi « Anan/m. Thus, one finds that ka \Z8nna 3 ps £}^ 2 . In other words, the de Broglie wavelength of 
the incoming atom must be shorter than the healing length of the cloud. For the sake of simplicity, we also 

impose an additional (less critical) requirement ka ^ 1 insuring that the s-wave scattering approximation is sufficient. 
Combining these two requirements, one obtains the condition 

C 1/2 <fca<l, (1) 

which has a wide range of validity as long as the gas parameter £ 1. 
One can represent the interaction Hamiltonian in the traditional form 

Hint = ^ / dx*t(x)*t( x )tf( x )tf( x ) , U = — , (2) 
2 J m 

Under the condition ([!]) the fast atoms behave almost like free particles, and their interaction with the rest of the gas 
can be treated perturbatively. Hence, it is reasonable to subdivide the total field \t into the low- and the high-energy 
parts, ip and ip', respectively: 

* = v+v', V' = 5> ke ' kx ' ( 3 ) 

k 

where the incident and the scattering states are described in terms of plane waves normalized to unit volume; ak 
destroys the high-energy particle with the momentum k. Substitution of Eq. (||) into Eq. (|J) and selection of the 
terms that describe the direct and the reversed processes yields 

H U = u o a ki a k 2 ak ^q + H ' c - ' ( 4 ) 

ki,k2,k 

where the Fourier component tp q = J dx exp(— iqx)^(x) of the field ^>(x) is introduced, and q = ki+k2 — k should 
be interpreted as the momentum transferred from (to) the cloud. 

First, let us consider the direct process. Then, in the initial state there are a fast atom with the momentum k 
and the cloud of N atoms. In the final state there are two fast atoms with the momenta ki, k2 and the cloud of 
N — 1 atoms. It is assumed that two outcoming atoms can be identified E3), so that the momentum transfer q 
and the energy transfer ui = (k\ + k 2 — k 2 )/2m can both be measured. Then, the corresponding double-differential 
cross-section for the direct process in the lowest order with respect to H[ nt is given by the Golden rule formula 



W(q , w) = Hip! J ^s{u-w fi ) J I d Xl dx 2 J°° dte^-**)-** t p(x 1 ,i;x 2 ,0) 



(5) 



Here p(xi, t\\ x 2 , t 2 ) = (0" t (xi,ti)^(x 2 ,t2)) [Note that unlike the uniformity in time, p(xi, ti', x 2 , t 2 ) = p(xi,ii — 
£ 2 ;x 2 ,0), the space uniformity cannot be, in general, assumed as long as a trapping potential exists]; ujfi = [q 2 /2 + 
qk + k\ — (q + k)ki]/m is the difference of the kinetic energies of the fast atoms in the final and initial states where 
the momentum of the second atom, k 2 , is fixed by the relation k 2 = q + k — ki. 

The integration over ki in Eq.(||) can be carried out explicitly. However, first we notice that the requirement of large 
k means that the values of q and ui, which are effectively selected by the correlator p(xi,i;x 2 ,0) in Eq. (H), satisfy 
the conditions q <C k and \u>\ <C k 2 /m. This immediately leads to the approximation S(oj — uJfi) ~ rn5(k{ — kki), 
which yields 

W(q,co) =4a 2 J J dxidx 2 J dte 1 ^-^- 1 ^ p(x u t;x 2 ,0) . (6) 
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We thus see that W(q, 10) is directly related to the dynamic correlator p(xi, t; x 2 , 0) which, contains rather rich 
information about the system, including, for one thing, the elementary excitation spectrum. Confining ourselves to 
the static correlations, described by OPDM p(xi,x 2 ) = p(xx,0;x2,0), we arrive at the even simpler relation 

W(q) = 87m 2 / / dx 1 dx 2 e lq ( Xl ^ X2 V(xi,x 2 ) (7) 



in terms of the differential cross-section W(q) = J du> W(q,oj). We note that such a simple correspondence between 
the OPDM and the crossection is due to the scattering length approximation of the two-particle interaction (see 
Eq.(|])). If the momenta of the fast atoms are larger than 1/a, no simple relation can be obtained in general, and 
only an estimate of the mean correlation length can be extracted from W(q). 

Now let us consider the time-reversed process. In this case the initial state consists of two particles with momenta 
ki and k 2 and N atoms in the cloud. In the final state there are one fast atom with momentum k and N + 1 atoms 
in the cloud. We note that the process under consideration can be viewed as the 4-wave mixing of the matter waves 
(see [ fL5| and references therein). In other words, two incoming beams of fast and coherent atoms aimed at the cloud 
may produce a third (outcoming) atomic beam. Its intensity J(n) per unit solid angle in the direction of some unit 
vector n is given by the Golden rule formula where the perturbation is the second term in Eq. (^) , and the finite state 
summation should be performed over the absolute value of k = nk. Finally, we arrived at the expression 

J(n) = 8g2 ^ nk2 Jduok{uj) J J dx^J 00 ^e lq ( xl - X2 )- Mt p W (xi,i;x 2 ,0) , (8) 

where the notation p^(xi, t; x 2 , 0) = (ip(x.i, £)V^(x 2 , 0)) is introduced; nk 17 nk 2 stand for the atomic densities of 
the incoming atomic beams; the momentum k = nfe(cj) of the outcoming atom carrying the energy excess 10 and the 
momentum excess q (if compared with the total energy and momentum carried by the two incoming atoms) is given 
by 



nfc(w) =q + ki +k 2 , k(u) = J k\ + k\ + 2mu>. (9) 

Note that when r c is large, q can be set to zero in Eq.(|9|). This implies, that w can be effectively selected by choosing 
values of ki, k 2 so that kik 2 = mu to satisfy the requirement of the energy and momentum conservation. 

Consider the structure of VK(q) and J(n) in the most characteristic cases. In the case of a pure Bose-Einstein 
condensate characterized by the condensate wave-function exp(— z/ii)<I>(x), one finds p(xi,x 2 ) = $*(xi)$(x 2 ) and 
p( r )(xi, t; x 2 , 0) = exp(— i/xi)p(x 2 , xi) . This yields 

W(q) = 87ra 2 |$ q | 2 , J(n) = 2n kl n k2 ™W(q), (10) 

irm 

where $ q = J dxe~ jqx $(x). In a trap having a center of symmetry, $(x) is characterized by a definite parity. Thus, 
$ q can be taken real. In this case ^(x) can be obtained from VF(q) ( measured either in the direct or the reversed 
process). 

An instructive case is the axially symmetric quantum vortex. The presence of a vortex carrying the vorticity 
I = ±1, 2, 3... in the center of the axially symmetric condensate drastically changes the scattering pattern. Indeed, in 
this situation, $ = exp(ilff)^J n(r, z), where 9 is the axial angle and n(r, z) stands for the axially symmetric condensate 
density as a function of the distances r, z perpendicular to the axis and along the axis, respectively. Accordingly, 



Eq. (10) shows that W(q) = for q directed along the vortex line. The differential cross-section becomes finite as 
long as there is a component qj^ of q perpendicular to the axis (or the vortex displaces from the condensate center) . 
In this case, W(q) ~ g 2Z for q± — > 0. 

Quite similar to the quantum vortex is the case of the supercurrent state of toroidal Bose condensate. The suppres- 
sion of W^(q) for q perpendicular to the plane of the torus permits the distinguishing of the supercurrent state from 
the currentlcss genuine ground state. 

Now let us consider the possibility of measuring the relative phase between two condensates by the discussed method. 
The problem of detecting a relative phase has been addressed many times in different contexts (see Refs. and 
references therein). The Bose field inside a split trap can now be represented as Jl7[ 

ip = oiV»i(x) + a 2 i/>i( x - Xq) (11) 

where a\ and a 2 destroy an atom in the condensate number 1 and 2, respectively; V'i(x) stands for the single-particle 
state of the first condensate; the single-particle state of the second condensate -01 (x — x ) is displaced by x from 
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the position of the first one. We assume that both states have no spatial overlap. Making use of the ansatz (|Tl]) in 
Eq.(^), one finds 

W{q) = 87ro 2 |Viq| 2 [<atai) + {a\a 2 ) + ((4a 2 ) e -* qx ° + h.c.)] (12) 

where ipiq stands for the Fourier transform of ipi (x). The means (a\ai) — Ni and {a\a 2 ) = N 2 represent the 
mean numbers N\, N 2 of atoms in the condensates. The interference effect is described by the crossterm (a} 122). If 
these two condensates are in a coherent superposition, a well defined relative phase tp can be introduced, so that 
(a\a 2 ) = exp(itp)y/NiN 2 . Accordingly, Eq. (|io|) acquires a form 

W(q)=2W 1 (q)(l+cos(^-qx )) ! (13) 

in the case Ni = N 2 = N. Here Wi(q) = 87ra 2 |-0i q | 2 A stands for the distribution which would be produced if only 
one condensate were present. Eq.([l3]) implies an extremely sharp dependence of VF(q) and of J(n) in Eq.([l0|) on ip. 

In the situation when the two condensates were initially in the number states lA^A^), no relative phase existed 
before the measurement. This implies that initially (a\a 2 ) = in Eq.(|l2]). Then, in the course of the scattering 
and measuring, the relative phase should be built in accordance with the general understanding p6|Jl7| . In order to 
describe the dynamics of the phase formation, the evolution of the correlator (a\a 2 ) should be found, where (...) is 
understood as the conditional quantum mechanical mean under the condition that a certain amount of the scattering 
events has been detected. 

In the absence of a genuine condensate (or for the above-the-condensate part of OPDM) one normally deals with 
the so-called quasi-homogeneous regime, when the off-diagonal correlation radius r c is much less then the scale of 
density variation given by, e.g., the size R c of the atomic cloud. In this case, it is reasonable to introduce the variables 
r = X2 — xi and R = (xi + x 2 )/2, and to represent VF(q) as 

W(q) = 8ira 2 J rfRp q (R) , p q (R) = J dr e^ qr p(r, R) . (14) 

In the quasi- homogeneous regime, the OPDM in the Wigner representation (p q (R)) has the semiclassical meaning 
of local (at the point R) distribution of the particle momentum q. Such interpretation is justified by the condition 
r c <C R c and by the definition of r c as a typical distance over which p(r, R) decays as a function of r. We thus see 
that in this case W(q) yields spatially averaged momentum distribution. Moreover, without contradiction with the 
uncertainty principle, this averaging can be partially (totally) removed by collimating the incident (both the incident 
and outcoming) beams. Even without removing the averaging, VF(q) contains valuable information about long-range 
correlations in the system, since the averaging does not affect the order-of-magnitude value of the correlation radius 
r c . If the system contains both the condensate and the quasi-homogeneous above-the-condensate fraction, the cross- 
section VF(q) is given by a combination of Eq.(|l0|) and Eq.([b|). 

Let us discuss general requirements for the atomic detector and the atomic sources. Currently, the only realistic 
detecting scheme we envision is that suggested in Ref. This method relics on imaging the atomic interference 

fringes by resonant light H. In our case, a typical angular distance between the fringes is given by a as l/r c k -C 1 
as long as the conditions (|J) and r c n -1 / 3 are fulfilled. Choosing n = 10 15 cm -3 , a — 3 • 10 _7 cm, one finds the 
gas parameter £ w 3 • 10 -5 . Then, in accordance with Eq. ([!]), k w 2 ■ 10 4 — 3 • 10 6 cm _1 and for r c w 10~ 3 cm the 
angular resolution should be a w 5 • 10~ 2 — 5 • 10 -4 . The presence of the trapping field characterized by the change 
U of the potential energy will distort the fringes by producing the uncertainty 8k ~ m\U\/k of the momentum. This 
introduces the limitation 8k < l/r c or \U\/fi < ka 2 /Jfr c ). For the chosen parameters, it implies \U\/fi < 1. A source 
of fast |l9f| and coherent atoms - the atomic laser |20[ - must have an uncertainty Ak of the incoming momenta less 
than l/r c . Thus, Ak/k = a. Recently this ratio has been achieved as small as 0.09 pl[ . We note, however, that this 
is still not enough for resolving r c comparable with typical condensate sizes. Nevertheless, detecting the initial stages 
of the ODLRO formation appears to be quite feasible because r c in this case can be as small as ~ 10~ 5 cm, which 
relieves the restriction on Ak/k as well as that on \U\/ fi by a factor of 10 2 . 

In conclusion, we have suggested a method of scattering fast atoms in a pure enough momentum state off a trapped 
atomic cloud in order to test directly the one-particle density matrix of this cloud. The differential cross-section of the 
inelastic process, when one incoming fast atom produces two fast ones, as well as of the reversed one allow measuring 
the correlation length of the local off-diagonal order. This gives, in particular, a powerful tool for testing different 
scenarios for the formation of the off-diagonal long-range order in the traps. This method can also be employed for 
detecting quantum vortices and supercurrent states, as well as the effect of quantum depletion of the condensate. A 
relative phase of the split-condensate can be measured as well. 
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